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ABSTRACT 

We investigate a scenario in which feedback from black-hole X-ray binaries (BHXBs) sometimes begins in¬ 
side young star clusters before strong supernova feedback. Those BHXBs could reduce the gas fraction inside 
embedded young clusters whilst maintaining virial equilibrium, which may help globular clusters (GCs) to stay 
bound when supernova-driven gas ejection subsequently occurs. Adopting a simple toy model with parame¬ 
ters guided by BHXB population models, we produce GC formation efficiencies consistent with empirically- 
inferred values. The metallicity dependence of BHXB formation could naturally explain why GC formation 
efficiency is higher at lower metallicity. For reasonable assumptions about that metallicity dependence, our 
toy model can produce a GC metallicity bimodality in some galaxies without a bimodality in the field-star 
metallicity distribution. 

Subject headings: binaries: close — globular clusters: general — X-rays: binaries 


1. INTRODUCTION 

The formation of bound Globular Clusters (GCs) takes 
place in an exceptionally complicated environment. A com¬ 
plete explanation for observed GC pop ulations likely involves 
multiple physical processes (se e, e.g., [Portegies Zwart et al.| 
|2010 [Longmore et al.||2014l, probably including galaxy 
mergei histories (see, e.g.,|Kruijssen|2014|). However, a com¬ 
mon simple picture is that proto-clusters become unbound 
when supernova feedback (SNF) ejects the remaining cluster 
gas, if the gas fraction at that time is too high (^ 0.7). If true, 
that might suggest that that GC formation is linked to regions 
of very high star-formation efficiency (^ 0.3). 

It has been argued that young clusters lose th eir gas too 
early to be consistent with that picture (see, e.g., |Longmore| 
|et al.||20T4l and references therein). However, the inferred 
ages of known gas-free clusters which are most similar to 
young GCs are still broadly consistent with the model we 
present ( |Basfian et al.||201?| l, especially given the likely un¬ 
certainty liTHetermmmgTEeir ages (see, e.g., [Schneider et al.| 
2014|l. Nonetheless, if proto-GCs are shown to always lose 


their gas before SNF occurs then our model would not work. 

In this Letter we examine whether X-ray binary (XB) feed¬ 
back (XBF) might affect the survival probabilites of young 
star clusters, since black-hole (BH) X Bs should sometimes 
form before significant SNF occurs (see|Justham & Schawin-| 
|ski||20I2| l. BHXBs may decrease the gas fraction inside the 
proto-cluster, and so could increase the effective star forma¬ 
tion efficiency before SNe eject the remaining gas. 

XBF has already been suggested to be especially impor- 


tant in the epoch of galaxy formation (see, e.g., [Glover & 

Brand 

2003||Power et al.|2009 

Cantalupo 2010[ 

Mirabel et al. 

"MTT 

Justham & Schawinski 

2012| Artale et a 

t.||2015>. Not 

only might the radiative feec 

back from XBs 

3e significant 

(Cantalupo [2010 Jeon et al. 

2012| Power et al. 201311, so 


might their kinetic output ( [Fender et al.|2005) [Heinz & Grimm| 
|2005[ l. Some X Bs directly drive very energetic kinetic out¬ 
flows (see, e.g., [Gallo et al.[[2005[ [Pakull ^ al.[[2010[ [Soria[ 
et al.[[2014|l. Stochastic galaxy-to-galaxy variation in XBF 
may help to explain some of the diversi ty in dwarf galaxy 
populations ( Justham & Schawinski|20f^ . 

Here we argue that the infiuence of BHXBs on the sur¬ 
vival of proto-GCs might explain why the GC formation ef- 
ficency - tjqq, i.e., the fraction of stellar mass which stays 
in bound clusters - apparently increases at lower m etallicity 
(see, e.g., Harris & Harris[[2()02[ [Forte et aL]|2005 i. Such a 
change in rjcc is necessary to explain why the GCs in mas- 
e, metal- rich galaxies are predominantly metal-poor ([Peng 


[et al.|2008[ ). 
resents both 1 


The present-day observed rycc presumably rep- 
formation-e poch processes and subse quent clus¬ 
ter disruption (see, e.g., [Gnedin & Ostriker| 1997] and Figure 
[T]|. Unless the other processes which affect ?7 gc strongly 
metallicity dependent then the model we present should be 
able to accommodate them. 

We do not e xamine the formati on of multiple stellar popu¬ 
lations in GCs. [Leigh et al.[([2013|l have done so when consid¬ 
ering single stellar-mass BHs which accrete from the proto¬ 
cluster gas. If BHXBs can affect the formation epoch of GCs, 
then they might somehow help to resolve this problemj^ 
Section|^presents the background to this scenario. Sections 
[^and|^investigate the idea using toy models. 


2. FEEDBACK FROM X-RAY BINARIES IN EMBEDDED YOUNG 
CLUSTERS 

A key property of BHXBs for proto-GC survival is that 
BHXBs can be formed before SNe eject the cluster gas. This 


' By contrast, as we were finalising this manuscript, [Cabrera-Ziri et al.[ 
(20151 speculated that XBs might have ejected gas from some young clusters, 
thereby reducing the opportunity for second-generation star formation. 
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^formation X (1 “ ^disruption) X ^survival — ^GC 


Star formation 
produces young, 
gas-rich star 
clusters. 



In some young clusters, 
feedback from BH XBs 
occurs for ~1 Myr before 
strong SN feedback does. 




The decreased gas fraction 
should help such clusters 
to stay bound. 



Clusters may be 
disrupted later 

-► 

by, e.g., tides. 


Overall process 
gives apparent 
present-day GC 
formation 
efficiency. 



At low 
metallicity, 
rjGc seems to 
increase. 


Figure 1. The overall formation efficiency of present-day GCs ( 1700 ) might be determined by multiple processes. If r)GC is Z-dependent, as observations 
indicate, this may mostly be contained within disruption, which we suggest XBF may explain. In addition, or alternatively, formation could be Z-dependent, 
perhaps indirectly (if, e.g., low-Z star formation tends to happen at higher intensity or density). Strong Z-dependence of Survival seems harder; it might rely on 
the different locations of the Z-rich and Z-poor cluster populations, since their ages are similar. 


is possible if BHs can form relatively quietly, from stars more 
mas sive than those which dominate the energetic S NF (see, 
e.g., Heger et al.|2003 Justham & Schawinski|2012 and refs, 
th erein). 

Rappaport et al.|p005[) find young BHXBs turn on «4 Myr 


after a starburst; their models match the observed luminous 
XB distribution in the Cartwheel galaxy. The expected time 
available for such BHXBs to act before strong SNF depends 
on whether BH formation “by fallback” is sufficiently quiet to 
avoid sudden gas ejection, or only formation of BHs by direct 
collapse. Those options roughly correspond to initial single¬ 
star masses above « 25Mq and « 4OM0, respectively (a more 
complete picture w ould be metallicity-dependent - see, e.g., 
Heger et al. 2003 1 - and binarity adds further complexity). 


Nonetheless, by the time SNe explosively eject gas from the 
proto-GC, one or more BHXBs might well have had wl-S 
Myr to reduce the gas fraction inside the cluster, which is a 
signficant qualitative difference between those proto-clusters 
in which an early BHXB does form an d those in which one 
does not Pustham & Schawinski|2012| l. 

Ther e are exceptions to the abov e broad statements. For ex¬ 
ample, iBelczyns^^&^^a^ (|200^ predict that ~0.1 per cent 
of NSs form within 5 Myr of st^ formation due to binary 
interactions, presumably in energetic SNe. Even rarer are 
BH-forming events probably associated with long-dura tion 
gamma-ray bursts (see, e.g., Podsiadlowski et al.|2004 and 
references therei n), and classes o f “super-luminous” super- 
novae (see, e.g., |GaTYam]|2012) l. Our toy models neglect 
these unusual events, although they may be sufficiently com¬ 
mon to affect detailed cluster demographics. 

Another complication is the possibility of age spreads larger 
than «1 Myr within cluster stellar populations. Since the 
expectat ion is that denser clusters should have smaller age 
spreads ( |Elmegreen]|2000[ |Tan et aL]|2006[ |Parmentier et al.| 
|2014| l, this may make gas ejection by BHXBs relatively less 
important for less dense clusters than for proto-GCs. 


2.1. Energetics and binding energy 

Eor a BHXB to significantly alter the distribution of gas 
within a protocluster, the energy input must be at least com¬ 
parable to the binding energy of the gas. Adopting a bind¬ 


ing energy Ebind for that gas of order GMMg^Ls/R=GM^ R, 

(where /gas = Mgas/M is the fraction of the total cluster mass 
- M - in gas, and R is a characteristic radius for the cluster) 
we estimate the binding energy of the gas as: 


Ebind 


10^° erg 


M \V/gas\ ( R y’ 
105 M^y 'v 1 ) viopcy 


( 1 ) 


Using those characteristic values for M, R and /gas leads to a 
ratio between Ebind Exbf - the energy input from XBE of 
luminosity Lxbf over time fxBF - of order: 


f^bind 1 f f-XBF \ f ^XBF ^ 

which suggests that a single luminous XB could significantly 
affect the gas within a young cluster, even for feedback effi¬ 
ciencies as low as a few percent. 


2.2. Column Density 

Another argument in favour of the importance of XBE in 
young GCs is that the gas-rich environments of embedded 
clusters might be Compton-thick. This would help to trap the 
accretion luminosity. We take a column density of lO^'^ cm“^ 
as sufficient for Compton-thickness, or a surface density Eh 
of slightly above 1 g cm“^. Scaling to this value, we hnd: 


1 g cm“2 ) \ IO^Mq J \ 10 pc / 


( 3 ) 


where Mh is the hydrogen gas mass. So for gas fractions of 
« 90% (99%), proto-clusters with initial stellar masses above 
IO^Mq (10" ^ Mq) m ay well be Compton-thick. 

Zezas et al. P002|l reported six heavily obscured luminous 


X-ray sources consistent with being inside dense, gas-rich, 
star-forming clouds in the Antennae. Those may well repre¬ 
sent less extreme examples of such embedded systems. 


2.3. BHXB formation probabilities and metallicity 
dependence 
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For studying GC survival, we only consider the most 
rapidly-formed BHXBs - those which switch on before strong 
SNR Here we assume that the formation of those pre-SNF 
BHXBs can be treated as a Poisson process. 

At solar met allicity we take the mean population from the 
calculations of Rappaport et al. (20051, using the predicted 


number of BHXBs more luminous than 10^® erg s ‘ at an age 
of «5 Myr. At that time, for a population which will even¬ 
tually produce 10^ core-collapse SNe, their least optimistic 
model predicts «1 such BHXB, and their more optimistic 
models «10 of them. We adopt an intermediate value of 3 
(simply normalised to IOOM 0 of stars per core-collapse SN), 
so models allow at least half an order-of-magnitude of free¬ 
dom in either direction in this BHXB frequency. 

At “low metallicity” - roughly below 0.1 Zq - we assume 
that the incidence of suitable luminous young BHXBs is 10 
times higher than at Zq. Observational evidence of a sim- 
ilar increase has become increasingly convincing (see, e.g.. 


Mapelli et al. 

2010 

IPrestwich et aL|2013||Brorby et al.|2014 

Douna et al. 


1 , and an increase was strongly expected 

by models (e.g., Belczynski et al.||2004 

Dray||20061 Linden 

et al.|2010 Justham & Schawinski|2012 

Fragos et al.|2013p. 


change in which single stars would form BHs, since other ef¬ 
fects are important in forming a system containing a BH in a 
close binary. In particular, early envelope loss - with respect 
to the nuclear evolution of the stellar core - can lead a massive 


star to leave a NS remnant rather than a BH (see, e.g.. Well- 

stein & Langer|1999[ Brown & Lee 20041 

Belczynski & Taam 

2008)1. However, at sufficiently low meta 

licity, massive stars 


bum helium before expanding to become giants, hence their 
common-envelope phases tend to happen later in their nuclear 
evolution, hence they are more likely to form BHs in close 
binaries. This effect may produce a sudden increase in suit¬ 
able BHXB formation below the appropriate threshold metal¬ 
licity (whi ch might be significan t for Q. Obs ervational re¬ 
sults fr om |Prestwich et al.| ( |2013| l and especially [Douna et al.| 
(|2015|l also suggest a sharp increase in the frequency below 
a transition metallicity (below 12-l-log(0/H) « 8 for the cal¬ 
ibrations used, though conversion to an absolute metallicity 
scale is non-trivial). 




Initial stellar mass: log(M/MQ) 

Figure 2. The upper panel shows the fraction of clusters which remain 
bound in our toy model (dash-dotted curves represent a ~ Zq population, 
solid curves represent the low-metallicity cluster population). The lower 
panel shows how those survival efficiencies modify a simple power-law 
ICMF, which is shown using a dashed line. Using light grey lines, the lower 
panel also illustrates the outcome if the ICMF becomes flat below 10“^ Mq of 
stars. 


3. SURVIVAL PROBABILITIES AND FORMATION EFFICIENCIES 

We now use toy models to demonstrate potential conse¬ 
quences of this scenario. These adopt a very simple set of 
assumptions, including that star formation efficiency is so low 
that SN-driven ejection of the remaining gas would normally 
unbind a young cluster. So a cluster survives to become a 
GC if and only if a BHXB acts hrst. In the terminology of 
Fig.[n we assume that all stars form in clusters (^formation = T 
thou^ some clusters may have very low mass), and that if 
a BHXB is formed before SNF the cluster remains bound 
(??XBF = (1 -??disruption))- 

3.1. Toy model without later evolution 

First we neglect later disruption processes (pxBF = Vcc)- 
Then then probability of a cluster remaining bound is the 
Poisson probability of that mass o f sta rs forming one or more 
BHXBs before SNF (taken from §2.3[ l. In the statistical limit, 
this probability is also the fraction of clusters which remain 
bound with that particular initial mass (and metallicity). 

Figure presents outcomes from such probability calcu¬ 
lations. Tne upper panel shows survival fractions for given 


initial cluster masses, and the lower panel the effect on a pop¬ 
ulation of clusters with a given initial cluster mass function 
(ICMF). The “initial cluster masses” are stellar masses, not 
including gas. Since the survival fraction is very low below 
the mass range which typically produces GCs we conclude 
that, if this toy model is inappropriate below the mass range 
which produces GCs, it would only moderately affect broad 
predictions. 

Figure compares integrated s urvival efficienci es to 
present-day inferred rjcc values (from [Forte et al.|2005| l. The 
upper and lower panels demonstrate the effect of varying the 
uncertain upper and lower bounds of the integral over the 
ICMF. The values predicted for rjQc, along with the difference 
between pcc for metal-rich and metal-poor clusters, appear 
at least order-of magnitude consistent with observationally- 
derived values in most of the parameter space. 

Where the predicted tjqq appears too high, this allows free¬ 
dom for unclustered star formation and additional GC disrup¬ 
tion mechanisms. So it seems preferable if the ICMF flattens 
or truncates at a few hundred solar masses or more. This con¬ 
straint would relax slightly if suitable BHXBs are more com- 
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3 4 5 6 7 8 

Initial cluster stellar mass: log(M/MQ) 



2.0 2.5 3.0 3.5 4.0 

Minimum relevant initial mass: log(M/M0) 


Figure 3. The upper panel shows the cumulative value of r?xBF for a toy 
model when integrating over an ICMF with a lower-limit of IO^Mq. The 
lower panel demonstrates how the overall, integrated 77 XBF is affected by the 
lower mass limit of a power-law ICMF. Both panels show curves for ~ z© 
and “low-Z” populations, as indicated. Solid curves assume the ICMF has 
a power-law slope of —2, as in Fig. broken curves adopt an ICMF slope 
of—2.5. Horizontal col oured lines represe nt inferred rjQc values for red and 
blue GCs in NGC 1399 JForte et al.|2005). The abscissa in the lower panel is 
not necessarily the minimum mass of the ICMF; it may approximate where 
the ICMF flattens from a power-law. 

mon than we have assumed; we repeat that in BHXB models 
there is at least half an order-of-magnitude of freedom to al¬ 
low this. 


3.2. Toy model including later evolution 

To our earlier assumptions we now add a parameterized ver- 
sion of later cluster evolution (due to|Fall & Zhang|2001| via 
Jordan et al.|2007 1. For this we use a Monte Carlo approach, 


randomly drawing clusters from a Schechter-function ICMF 
and deciding whether each cluster remains bound based on 
the same BHXB formation probabilities as above. Each of 
the surviving clusters is then assumed t o lose a mass - de¬ 
noted by A - after the formation epoch ( Fall & ZhangpOOT 
[Jordan et al.|2007) l. 

Due to the large dynamic ranges, for the distributions in 
Figure]^ we drew 10" primordial clusters from the ICMF. 
Over the mass range in Fig. assuming A = lO^M© reduces 
rjcc by a factor of a few. This reproduces the GCMF turnover 
within our earlier model uncertainties. 



Stellar mass; log(M/Mg) 


Figure 4. The black solid histogram represents a Schechter-function ICMF, 
with an exponential cutoff scale of IO^Mq. The dotted curves (labelled 
A = 0) represent model predictions immediately after gas ejection. The solid 
coloured histograms present the mass functions after later mass loss (A, as 
labelled). The dashed purple lines show inferred slope s for metal-rich and 
metal-poor cluster populations from|Jordan et al.H2007). 


3.3. The upper slope of the cluster luminosity function 

Jordan et al.| (|2007)l found that the bright end of the GC 
luminosity function is steeper for GC populations around less 
massive host galaxies (which correlates with more metal-poor 
GC populati ons). Figure |4| displa ys the power law index (/3) 
inferred by [Jordan et al'| ( |2007] l for the upper mass func¬ 
tion of those presumed metal-rich and metal-poor populations 
(j3 = 1.7 and 3.0, respectively). The correspondence with this 
model seems at least intriguing. 

This qualitative effect also arises in our simplest toy model, 
as Fig. [^predicts that the upper mass function of bound GCs 
should^ steeper for low-metallicity clusters. 


4. CLUSTER METALLICITY DISTRIBUTIONS 

If the formation of appropriate BHXBs has a suitable metal- 
licity dependence, then this scenario might explain the metal- 
licity bimodality of Galactic GCs. It is not certain whether the 
commonly-observed colour bimodality of extragalactic GCs 
indicates a similar m etallicity bimodality ( [Yoon et ar][2006[ 
Blakeslee et al.|2012|l. A metallicity-dependent rjQciZ) would 
produce a GC metallicity bimodality when the host galaxy’s 
stellar metallicity distribution p,(Z) was appropriate - i.e. only 
if /r(Z) X rjGciZ) leads to a bimodal distribution. Unfortu¬ 
nately, robust niZ) distributions for entire galaxies, including 
their halo populations, are difficult to obtain. 

Here we show that a simple rjQciZ) function, guided by our 
model, could reproduce a GC metallicity bimodality from a 
non-bimodal ^(Z). We assume two regimes of constant rycc 
- one at low Z and one around Zq - matched by a transition 
region, expressed as: 




T-log[f- 


^o/ J 


(4) 


where R represents the amount by which log(r]Gc) increases 
at low-Z, S gives the “sharpness” of the transition between the 
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Figure 5. For given star formation Z-distributions (/^(Z); the dotted curves), we demonstrate how a toy model may recover a bimodal distribution in cluster 
metallicity (solid curves). Upper row: The fixed p{Z) is intended to be broadly representative of common p{Z) distributions. Each panel varies one of the 
model parameters, with different curves labelled with the relevant par ameter value. The dashed curve in the center panel shows the outcome when applying the 
metallicity dependence for generic HMXBs from|Fragos et al.H2013J. Lower row: For fixed parameters {R = 0.7, 5 = 5, T = 1), we show results from varying 
p{Z)\ in the left-hand panel, the black curves are tor the same /i(Z) as in the top row. 


z© and low-Z regimes, and T indicates the metallicity around 
which that transition occurs. 

For the upper row of Figure 1^ we assume a con¬ 
structed to a pproximate the /r(Z ) distribution of NGC 1399 
in figure 5 of Forte et ah] (2005 1 . There we keep p,(Z) hxed, 
and show cluster metallicity distributions for different values 
of R, S and T. Our adopted values of R and T are similar to 
those which we have previously sugg ested are expected (i.e., 
« 1 and « -1, respectively; see ^2.3[ ). From this /i(Z) we can 
recover a bimodality in cluster metallicity, and the qualitative 
result is not affected by sm all variations i n S. The gradual 
metallicity dependence from |Fragos et al.| ( |2013| l for generic 
high-mass XBs would not produce a bimodality, as shown in 
the central panel. However, the systems involved in our sce¬ 
nario are a very specific subset of young BHXBs. Moreover, 
there are in dic ations that a sharp t ransition is at least plausible 
(see § |2.3| and |Douna et al.|2015l l. 

The bottom row of Figure demonstrates changing /r(Z) 
with hxed model parameters (R = 0.7, 5 = 5, T = 1). This 
illustrates that the model does not predict that every galaxy 


must hav e a bimodal cluster m etallicity distribution (for which 
see, e.g., |Usher et al. 112012)1, and that the peaks of the GC 
metallicity distribution are affected by the underlying /r(Z). 

When future observations provide reliable p,(Z) distribu¬ 
tions for a signihcant sample of galaxies, this should enable 
strong constraints on our model. If Eq.|^ provides a good 
functional form for rjcciZ) - and if rjcc is only a function 
of metallicity - then the parameters R, S and T should be 
universal constants. However, we expect that reality is more 
complicated. For example, this simplihed model assumes that 
VxBF - and hence pcc - has only two regimes (Zq and low- 
Z). At higher metallicities pxBF might well decrease further. 
Indeed, for metallicities signihcantly above solar then sin- 
gle stars might no longer leave BH remnants at all (see, e.g., 
Heger et al. 2003 ] l. 


5. PROTON-CAPTURE REACTIONS AND ABUNDANCE 
ANTICORRELATIONS 

Since the Na-O and Mg-Al anticorrelations are exclusively 
associated with star clusters which stay bound, it would be 
elegant if these anticorrelations could be explained by the 
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BHXBs which - in this model - allow clusters to stay bound. 
This speculation arose since the relevant nuclear reactions are 


high-e nergy proton capture reactions (see, e.g., Izzard et al. 
|2007| l, and many BHX Bs may be prolific sources of high- 
energy protons (see, e.g.,|Heinz & Sunyaev|2002l|Fender et al.l 
1 ^ . - 


Unfortunately, even for assumptions which we expect are 
extremely optimistic, we have been unable to convince our¬ 
selves that a single BHXB would explain the observed abun¬ 
dance anomalies in an entire GC. 

Nonetheless, variations on this mechanism may deserve 
further examination, perhaps involving other proton accelera¬ 
tors in these dense stellar systems. 


6. CONCLUSIONS 

We have investigated a scenario in which young stellar clus¬ 
ters are able to survive as bound GCs only when feedback 
from BHXBs gradually decreases the gas fraction inside the 
cluster before energetic SNe suddenly eject the remainder of 
the gas. Several potential complications are ignored by our 
toy model. However, as long as the other processes which 
affect 77 gc are largely metallicity-independent, this scenario 
may well have sufficient freedom to include them without los¬ 
ing its positive features. A minimal set of assumptions, com¬ 
bined with input numbers estimated from population mod¬ 
els for suitable BHXBs, produces predictions consistent with 
observationally-inferred rycc^ including the increase in pcc at 
low metallicity. We therefore suggest that more detailed study 
of early feedback from BHXBs in proto-GCs may be impor¬ 
tant for trying to understand GC populations. 
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